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ABSTRACT 

We study the impact of early dark energy (EDE) cosmologies on galaxy properties by coupling 
high-resolution numerical simulations with semi-analytic modeling (SAM) of galaxy forma- 
tion and evolution. EDE models are characterized by a non-vanishing high-redshift contribu- 
tion of dark energy, producing an earlier growth of structures and a modification of large-scale 
structure evolution. They can be viewed as typical representatives of non-standard dark energy 
models in which only the expansion history is modified, and hence the impact on galaxy for- 
mation is indirect. We show that in EDE cosmologies the predicted space density of galaxies 
is enhanced at all scales with respect to the standard ACDM scenario, and the corresponding 
cosmic star formation history and stellar mass density is increased at high-redshift. We com- 
pare these results with a set of theoretical predictions obtained with alternative SAMs applied 
to our reference ACDM simulation, yielding a rough measure of the systematic uncertainty of 
the models. We find that the modifications in galaxy properties induced by EDE cosmologies 
are of the same order of magnitude as intra-SAM variations for a standard ACDM realization 
(unless rather extreme EDE models are considered), suggesting that is difficult to use such 
predictions alone to disentangle between different cosmological scenarios. However, when 
independent information on the underlying properties of host dark matter haloes is included, 
the SAM predictions on galaxy bias may provide important clues on the expansion history 
and the equation-of-state evolution. 

Key words: early Universe - cosmology: theory - cosmology; cosmological parameters - 
galaxies: formation - galaxies: evolution - galaxies: fundamental properties 
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1 INTRODUCTION 

The last decade has seen considerably advances in our understand- 
ing of the properties of our Universe as a whole, in particular thanks 
to the accurate meas urement of the most important cosmological 
parameters (see e.g. iKomatsu et al]|2009l and references herein). 
One of the most surprising discoveries of this "precision cosmol- 
ogy" epoch is that some unknown form of Dark Energy (DE, here- 
after) accounts for more than 70% of the energy density of the 
present-day Univers e, and is responsible fo r its accelerated expan- 
sion today (see e.g.. IPerlmutter et al]| 1 9991) . The physical nature of 
DE, together with its origin and time evolution, is one of the most 
enigmatic puzzles in modern cosmology. 

In the standard ACDM cosmological model, DE is treated as 
a classic cosmological constant (following Einstein's original con- 
jecture) where a homogeneous and static energy density fills the 
whole Universe. This simple assumption, however, gives rise to a 
number of theoretical problems, such as the high degree of "fine- 
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tuning" required to accommodate the present day value of Q,a (see 
I Weinberg! 19891 for a review). For this reason, many alternative sce- 
narios to explain DE and the accelerated expansion have been pro- 
posed, ranging from scalar field models such as quintessence to rad- 
ical modifications of the laws of gravity. Most observational efforts 
currently concentrate on constrai ning effective parametrisations of 
the D E equation of state (see e.g. lWetterichlll988l : lRatra & Peeblei 
Il988l) . which is adequate for 'non-coupled' dark energy models in 
which the dark energy can be approximated as uniform and influ- 
ences structure formation only through a modification of the cosmic 
expansion rate. 

An interesting sub-class of these models are scenarios that 
involve a non-negligible DE contribution at early times dur- 
ing recombination and primordial structure formation, which 
also implies non-ne gligible modificati ons of the cosmic mi- 
crowav e background i Doran et alj|200lh , of big-bang nucleosyn- 
thesis {MUU^etaL 2004 ) and of large-scale structure forma- 
tion teartelmann et alj|2006l) . The non-linear structure formation 
predicted in such early dark energy (EDE) models has been 
studied through high-resolution A-body simulation l lBaldi|[201^ : 
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iGrossi & SDrii^l2009l) . In particular, the impact on the statisti- 
cal properties of dark matter (DM) haloes (such as the abundance 
of high-redshift galaxy clusters) as a function of cosmic time has 
been compared with the corresponding predictions for the stan- 
dard ACDM cosmology, with the aim of identifying observational 
tests that are able to disentangle between the differe nt cosmologies 
based on future surveys (like the EUCLID satellite. iLaureiis et alj 
l20Tlh . 

In contrast, the influence of modified DE models on the prop- 
erties and the evolution of galaxy populations has not yet been 
explored in full detail, even though many cosmological tests for 
constraining the DE equation of state ultimately rely on a pre- 
cise understanding of how galaxies trace mass. This can in part 
be understood as a result of our limited present understanding of 
galaxy formation even in the ACDM cosmology, which is ham- 
pered by several long-standing discrepancies between the predic- 
tions of theoretical models and observations, as seen, for example, 
in t he redshift evolution of the galaxy stellar mass function (see 
e.g. lFontanot et alj20()9l). the properties of the Milky Way satellites 
(see e.g. Maccio et al. 2010l : iBovlan-Kolc hin et al."2012'). the low 
baryon fraction in galaxy clusters (see e.g. McCarthy et al. 20o3), 
or in the shallow DM profiles associate d with differen t galaxy pop- 
ulations ("cusp-core" problem, see e.g. lMoorelll994h . This in turn 
fuels some interest in exploring both the possibility that these dis- 
crepancies may be reduced by assuming an evolving DE, and in 
finding observational tests based on statistical galaxy properties 
that could potentially distinguish between different DE scenarios 
(which one may also hope to be easier to perform compared with 
tests working purely in the "dark sector"). 

The catch of course is that galaxy evolution is a complex pro- 
cess, involving a non-linear blend of many different physical pro- 
cesses acting on the baryonic gas. All models of galaxy formation 
on cosmological scales (both semi-analytic models, SAMs here- 
after, and numerical simulations alike) need to simplify this intrin- 
sic complexity by means of quite coarse, yet physically grounded, 
analytic approximations (in SAMs) or sub-grid models (in simu- 
lations). In SAMs, these analytic approximations are meant to de- 
scribe the physical processes acting on the baryonic gas (such as 
gas cooling, star formation and feedback) as a function of the phys- 
ical properties of model galaxies (like their cold gas and stellar con- 
tent). This then yields a convenient parametrisation of the physics, 
which is calibrated against a small subset of low-redshift observa- 
tions. The ability of modern models of this kind to successfully 
match a large and diverse set of observations not used in the cali- 
bration can be viewed as a powerful confirmation of the viability of 
the basic paradigm of hierarchical galaxy formation. 

However, it has also been shown that the SAM approach en- 
tails a sig nificant level of degene racies among the different param- 
eters (e.g. iHenriques et"Zll2009l) , and the theoretical uncertainty is 
exacerbated by the fact that different models tend to adopt different 
parametrizations for the physical processes acting on the baryonic 
gas. Even though the comparison of different model predictions 
shows a reassuring level of consistency in many cases (see e.g. 
iFontanot et al.ll2003) . there is hence a certain degree of systematic 
uncertainty in this approach. Therefore, in order to assess the role 
of modified DE on galaxy evolution and, viceversa, to use galaxy 
formation to constrain DE models, it is important to not only quan- 
tify the impact of modified cosmologies on galaxy formation within 
a specific model, but also to check how the size of the changes in 
the model predictions compares to the intra-model variance for a 
fixed ACDM cosmology. 

This paper is organized as follows. In Section[2l we introduce 
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the cosmological numerical simulations and semi-analytic models 
we use in our analysis. We then compare the various predictions in 
Section|3] Finally, we discuss our conclusions in Section|4] 



2 MODELS 

2.1 Dark energy parametrisation 

In this paper, we focus on a specific class of DE cosmologies, called 
Early Dark Energy (EDE) models. In these cosmologies, DE con- 
stitutes an observable fraction of the total energy density at the time 
of matter-radiation equality (in contrast to the cosmological con- 
stant scenario). A natural outcome of these models is to predict an 
earlier formation of structures with respect to the ACDM cosmol- 
ogy for an equal amplitude of the present-day clustering strength. In 
fact, in EDE models the high-z cluster population grows consider- 
ably relative to ACDM, helped also by a lowere d value of the crit- 
ical l inear density contrast needed for collapse teartelmann et al.l 
1 20061) : as a result, EDE models predict a slower evolution of the 
halo population than the standard ACDM cosmology. 

In particular , we consider similar EDE mod els as in 
iGrossi fe Sprin"^ ( |2009|) . as originally introduced by IWetterictil 



20041) . These models are characterized by a low but non- vanishing 



dark energy density at early times, and by a low redshift value of 
Q.A consistent with CMB constraints. The equation of state param- 
eter wi^z) varies with time according to the equation 

wo 



w{z) = 

where 

fe = 



3 wo 



(1) 



(2) 



In the previous equations, wq and Qdc.o = 1 — f2m,o represent 
the present-day equation of state parameter and amount of dark 
energy, respectively, while f2do,c gives the average energy density 
parameter at early times. For sufficiently low f2dc,o, EDE models 
reproduce the accelerated cosmic expansion observed in the local 
Universe. 

We examine a set of different EDE models by changing wo 
and f2dc.c (Table[Tll. In particular, we vary the f2dc,c parameter over 
two orders of magnitude; it is worth stressing, however, that in our 
EDE5 model we assum e a value for Qdc,c whic h is rather extreme 
for this class of models ( iHollenstein et al.ll2009l) . 

2.2 Numerical simulations 

We have performed a series of numerical A'^-body simulations of 
DM-only runs for all our EDE models. For comparison, we also 
run a simulation in a standard ACDM cosmology configuration and 
we consider a model with constant equation of state w = —0.6 
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Figure 1. Redshift evolution of the galaxy stellar mass function in different cosmological models, and for different SAMs in the same ACDM model. Various 
line st yles and colours distinguish the different models, as labeled in the legend. Grey points refer to the stellar mass function compilation from lFontanot et al] 
see references herein). 



(CWDE) as well. In all simulations, we assume a flat universe with 
matter density parameter f7m ~ 0.25, Hubble parameter h — 0.73, 
and Gaussian density fluctuations with a scale-invariant primordial 
power spectrum with spectral index n = 1. The normalization of 
the linearly extrapolated z — power spectrum is taken to be erg = 
0.9 for all simulations. 

We generate initial condition for all the cosmological models 

using the N-GENIC code and have run the simulations using the cos- 

I : — — 1 

mological code GADGET-3 (last described in'Springel 2005). For 
both codes we employ the versions modified bv .Grossi & Springeli 



( l2009l) to account for EDE cosmologies with a redshift-dependent 
equation of state. All simulations have been run using 432^ par- 
ticles in periodic boxes 100/i~^Mpc on a side, corresponding to 
a mass resolution of 1.18 x 10^ h~^MQ. For each box, 64 sim- 
ulation snapshots w ere saved, at the sam e redshifts used in the 
Millennium project ( ISpringeletal]|2005l) . facilitating a straight- 
forward comparison. Group catalogues have been constructed us- 
ing the friend-of-friend (EOF) algorithm with a linking length of 
0.2 in units of the mean particle separation. Each group has then 
been decomposed into gravitationally bound substructures using 
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Figure 2. Redshift evolution of (a) tlie stellar mass density, and (b) the cos- 
mic star formatio n rate density . Grey points refer to a stellar mass density 
compilation from lSantini et al.1 )2012l see references herein) and to a cos- 
mic star formation density compilation from lHop kins 1 2004;, see references 
herein), respectively. Our models are shown with the same lines types and 
colours as in Fig.[T] as labeled. 



SUBFIND l lSpringel et al.ll200ll) . The resulting subhaloes are then 
used to construct m erger history trees as explained in detail in 
ISpringel et al] l l2005l) . Only subhaloes that retain at least 20 bound 
particles after the gravitational unbinding procedure are kept for 
the tree construction, thus implying a subhalo detection threshold 
of 2.36 X lO^°/i"^Af0. 

As a consistency check, we have measured the matter power 
spectrum of the different simulation m odels at various redshifts , 
finding consistent results with respect to lGrossi & Springell l [2009l , 
their Fig. 5). Similarly, our results for the halo mass functions re- 
produce their findings. 



2.3 Semi-analytic models 

In this paper, we mainly co nsider the most recent implemen- 
tation of the Munich model jGuo et al] 1201 ih. based on th e L- 
GALAXIES code originally developed bv lSpringel et al. l( l2005l) .We 
also consider two p revio us versions of the same cod e, namely the 
ICroton et al] j2006l) and iDe Lucia & BlaizoJ ( |2007|) models. The 
use of different SAMs implemented on top of the same underly- 
ing code structure allows an assessment of the intra-model variance 
induced by different choices for the approximation of galaxy for- 
mation physics. This yields a course estimate of systematic model 
uncertainties, which can in turn be compared to the size of changes 
in the predicted galaxy properties due to EDE cosmologies. The 
chosen set of SAM models is particularly suitable for this ap- 
proach as it represents a coherent set of models designed to work 
on Millennium-like merger trees, such that complications such as 
merger tree conversions can be avoided. 

In the interest of conciseness, we refrain in the following from 




Figure 3. Deviation of the galax y auto-correlation functions in our model 
cosmologies with respect to the IGuo et afl j201 ih model applied to the 
ACDM simulation (5ga)- The lines types and colours are the same as in 
Fig.E 



discussing all the details of the modeling of different physical pro- 
cesses in each of the SAMs (we refer the reader to the original 
papers for a full specification). However, we want to point out 
where the main differences between the three models lie: (a) in 
the treatment of dynamical friction and merge r times, the stellar 
ini tial mass function and the dust model (from lCroton et al]|2006l 
to lOe Lucia & Blaizolll2007l) : (b) in the modeling of supernovae 
feedback, the treatment of s atellite galaxy evolution, ti da l stripping 
and m ergers (improved from lPe Lucia & Blaizoll2007l to lGuo et"aD 
I2OIII) . In all cases, these major changes in the physical recipes in- 
volved a re-calibration of the most relevant model para meters. 

W e note that we have modified the code used in the lOuo et al.l 
l|) model in order to have it run self-consistently on EDE cos- 
mologies, by including the EDE contribution in the expression for 
the Hubble expansion rate: 



H{a) = Ho^^ + fidco e-^/[l+'^("'l'''"^ 



(3) 



This equation is a generalisation of the usual Hubble function defi- 
nition employed in L-GALAXIES , and it naturally simplifies to the 
correct forms when dealing with the CWDE and ACDM cosmolo- 
gies iw{a) — —1). In the fo llowing, whenever we discuss results 
based on the lOuo et al 11) model, we actually refer to this mod- 
ified code. 

All SAMs have been calibrated by requiring them to repro- 
duce a well-defined set of low-redshift reference observations. In 
this paper, we do not discuss possible changes in the model cali- 
brations, rather we prefer to retain the original parameter choices, 
in order to allow a direct comparison to the published models and to 
highlight differences induced by changes in the cosmology alone, 
for physics models that are kept fixed. As far as the EDE and 
CWDE cosmologies are concerned, this implies that these mod- 
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Figure 4. Redshift evolution of the galaxy bias in different cosmological models. The lines types and colours are as in Fig.[T] 



els are not necessarily tuned to perform best, as in the ACDM case. 
Nonetheless, since we require our numerical simulations to provide 
consistent statistical properties for z — DM haloes by construc- 
tion, we expect the variations between SAM predictions in the dif- 
ferent cosmologies to be small at the lowest redshifts. I n any case, 
our ch oice of keeping the parameter choices used in the lGuo et al.l 
( I2OI ih model the same also in different dark energy cosmologies is 
instrumental for testing the effective change in model predictions 
due to a variation of the assumed Hubble function alone, a distinc- 
tive feature of EDE models. 



3 RESULTS & DISCUSSION 

We begin by comparing the redshift evolution of some of the most 
basic global properties of the galaxy populations predicted by our 
SAMs, namely the galaxy stellar mass function (Fig. Q}, and the 
evolution of the cosmic star formation rate and stellar mass density 
(Fig.O. 

In all the figures, the black lines show the different L- 
GALAXIES mode ls in a ACDM co smology. In particular, the solid 
lines refer to the lOuo et al] ( 1201 ih implemen tation of the model, 
while dash ed and dotted line s refer to the iDe Lucia & Blaizoj 
( I2OO7I) and lCroton et"ai] ( l2006i) versions, respectively. The differ- 
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ences between the predictions of thie different SAMs applied to 
the same cosmological simulations are mainly driven by the dif- 
ferent cali bration sets emplo yed by the different authors. For ex- 
ample, the IGuo et af] ( 1201 ll) model has been especially designed 
to reproduce the low-mass end of the stellar mass function at 
z = 0. Co l oured lines, on the other hand, give the predictions of the 
iGuo et al. I l201lh model in different cosmologies. Red lines refer to 
the CWDE simulation, whereas blue, green and yellow lines refer 
to the different EDE realizations, as labelled in the legend. 

We point out that the agreement of the lGuo et alj jlOllh pre- 
dictions among the different cosmologies is satisfactory at 2: = 0. 
When a comparison to observational constraints is made, the model 
predictions have been convolved with an adequate estimate of 
the typical observational error (e.g. a lognormal error distribution 
with amplitude 0.25 and 0.3 for stellar ma sses and star formation 
rates respectively, s ee lFontanot et al1l2009l) . The predictions of the 
ICroton et al] ( I2OO6I) model have been converted from a Salpeter to 
a Chabrier IMF by assuming a constant shift of 0.25 dex in stellar 
mass and 0.176 dex in star formation rate. 

As expected, the earlier structure formation epoch in EDE cos- 
mologies directly translates into an increase of the space density of 
galaxies at all redshifts (Fig.[T). The strength of this effect and the 
redshift range involved increase with the Sldc.o value, as then the 
dark energy contribution to the total energy density of the Universe 
becomes more relevant. Nonetheless, it is worth noting that this in- 
crease in space density is not limited to the massive end of the stel- 
lar mass function, but is present at all mass scales. The well-known 
problem of an overprediction of intermediat e-to-low mass galax- 
ies in theoretical models of galaxy formation jFontanot et al.l200^ 
is thus not reduced, but even enhanced in EDE cosmologies. Fol- 
lowing this argument, we also expect that a number of known ten- 
sions between model predictions and the properties of galaxies in 
the Milky-Way environment should not be solved but exacerbated, 
although dedicated higher resolution simulations are needed to firm 
up this conclusion. However, when we consider the predictions of 
the other two L-GALAXIES versions based on the same ACDM cos- 
mology, we find a small spread in the predictions of different SAM 
models computed for the same cosmology. This "systematic uncer- 
tainty" in the theoretical modelling is of the same order as the size 
of the differences seen in the predictions of a single SAM model 
applied to different cosmologies. 

Similar conclusions can be drawn by considering integrated 
quantities, like the cosmic stellar mass density and mean star for- 
mation rate density (Fig. |2l. In these quantities, the different red- 
shift evolution among different cosmologies and models is even 
more clear: out to redshift 2 ~ 4, it is difficult to separate the im- 
pact of different cosmologies from the effect of a different model- 
ing of physical processes and/or different model calibrations, while 
at higher redshift, EDE models effectively predict somewhat higher 
star formation rates and accumulated stellar mass density. Anyhow, 
also in this case, only rather extreme models differ significantly 
from the locus of considered ACDM SAMs. 

The galaxy auto-correlation function ^gai represents another 
potentially useful discriminant between different cosmologies. We 
consider auto-correlation functions corresponding to different in- 
tervals of stellar mass. The predicted galaxy auto-correlation func- 
tions for our models show good consistency among each other over 
a wide redshift range, so they provide no direct way to disentan- 
gle between different cosmologies. In Fig. [3] we show the resulting 
aut o-correlation func tions normalised to the predicted values for 
the IGuo et al] l l201lh model in the ACDM realization (^ga), for 
all galaxies more massive than IO^Mq and in two different mass 



bins. Similar trends are seen in all cosmologies, with no relevant 
deviations from the reference model in most cases. Most not ably, 
the strongest deviat i ons ar e seen for the ICroton et al.l ( l2006h and 
iDe Lucia & BlaizotI JioOTl) predictions at very small scales. This is 
cl early due to the d ifferent treatment of satellite galaxy evolution 
m IGuo et alj l [20Tlh . For the theoretical predictions on these small 
scales, the physical mechanisms acting on baryonic gas dominate 
over the cosmological effects. 

Nonetheless, these results can provide interesting insights 
when coupled with additional information on the underlying dis- 
tribution of DM particles. In particular, for each cosmology, we 
have computed the auto-correlation function ^dm of the mass for 
a randomly selected subsample of DM particles (corresponding 
to 1% of the total particles in the cosmological box, for compu- 
tational convenience - a larger fraction would not change the re- 
sults). As expected, the redshift evolution of ^dm shows large de- 
viations among different cosmologies. We combine these comple- 
mentary informations into an estimate of the galaxy bias, defined as 
Cgai/Cdm (Fig. 11). In all stellar mass intervals we consider, it is in 
principle possible to distinguish between the different cosmological 
models with the help of the bias, by considering the intermediate to 
large scales (i.e. larger than a few tenths of a Mpc) at 2 > 0. In 
fact, at these scales, there is a concordance of the predictions for 
the bias from the different ACDM based SAMs, while, at the same 
time, the deviations of EDE models from a standard ACDM cos- 
mology become progressively larger at increasing redshift and for 
larger stellar mass. At smaller scales, the intrinsic intra-model dis- 
persion between SAM predictions due to the different treatments 
of the satellite galaxy physics hampers firm cosmological conclu- 
sions. 



4 CONCLUSIONS 

In this paper we discuss the expected impact of Early Dark En- 
ergy (EDE) cosmologies on the properties of galaxy populations, as 
predicted by semi-analytic models. We cons ider the latest version 
of the L-GALAXIES model jGuo et al]|201lh . modified to account 
for time-dependent variations of the equation of state parameter w. 
We a l so consider earlier version s of the same SAM jCroton et al.l 
l2006l : lDe Lucia & Blaizoll2007l) to compare the size of differences 
in galactic properties induced by changes in the underlying dark en- 
ergy model with the intra-model variance due to differences in the 
modeling of the baryonic physics. The extension of our approach 
to consider other c osmological models with alternative dark energy 
scenarios (see e.g. lBaldill2012l) is in principle straightforward, pro- 
vided high-resolution N-body simulations of such scenarios can 
be calculated. The latter has very recently become possible even 
for more compli cated theories of gravity, such as DGP or ,f(R)- 
gravi t y (see e. g. 



ty, „ , 

Oyaizul I2OO8I : ISchmidlll2009l : iKhourv & WvmanI 
2|), such that our method can be fruitfully used 



it v (see 

l2009l : lLi et all^Oi; 
to connect such scenarios with observations of galaxies. We plan to 
expand our present results into this direction in forthcoming work. 

Our results highlight that EDE cosmologies lead to important 
modifications in the galaxy properties with respect to a standard 
ACDM universe. Nonetheless, they also show that these deviations 
are of the same order of magnitude as those induced by different as- 
sumptions in the modeling of physical properties and/or by parame- 
ter r ecalibrations in a homologous set of ACDM-based SAMs (see 
also IGuo et al. ||20l2). Stronger effects are seen at higher redshift 
(2 > 4), but these redshifts correspond to a cosmic epoch where di- 
rect observational constraints on galaxies are extremely difficult to 
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obtain and where we expect SAM predictions to be comparatively 
uncertain (given that their usual calibration samples are selected 
at lower redshifts). This behaviour is mainly due to our still lim- 
ited understanding of the physical properties acting in the baryonic 
sector, and to the high level of degeneracy between the analytic ap- 
proximations employed in different SAMs. We thus conclude that 
predicted galaxy properties alone provide only weak constraints on 
disentangling a standard ACDM cosmology from an Early Dark 
Energy model. These results are consistent with studies that anal- 
ysed the response of SAMs against variations of cosmological pa- 
rameters (coffesponding to the best constraints from different data 
releases o f the WMAP satellite) in the ACDM Universe (see e.g 
IWang et al.. 2008 : Quo et al., 201% . 

We have also shown that some of these degeneracies are bro- 
ken if additional information on the properties and distribution of 
the underlying DM field became available. The combination of the 
stronger evolutionary patterns expected for the "dark sector" and 
the definition of galactic populations with suitably stable proper- 
ties in SAMs indeed leads to the identification of reliable cosmo- 
logical tests. In particular, we show the potential of a measurement 
of galaxy bias on scales larger than a few tenths of Mpc as a cos- 
mological discriminant. Future observational efforts aimed at con- 
straining w and its redshift evolution, like the EUCLID mission 
jUaureiis et alj20l"lh . are indeed devised in order to provide, at the 
same time, constraints on the DM distribution (via weak lensing) 
and on the properties of large galaxy populations (via spectroscopy 
and photometry). We stress however, that this mission is designed 
to primarly focus on the z ~ 1 diagnostics, where the predicted ef- 
fects of EDE cosmologies on galaxy bias are smaller than at higher 
redshifts. 

While we confirm that the combination of these different 
probes is indeed one of the most promising methods to achieve 
constraints on the dark energy component of our Universe, we also 
would like to emphasise the power of the particular simulation ap- 
proach taken here, which is able to accurately predict galaxy bias 
in non-standard cosmologies as a function of scale, epoch, and 
galaxy type. The combination of high-resolution dark matter simu- 
lations and semi-analytic galaxy formation models used in this pa- 
per should also be highly useful to explore and understand generic 
effects of galaxy formation on cosmological constraints, for in- 
stance tho se derived from BAO or redshift space distortion mea- 
surements l lAngulol20l3) . This would ultimately lead to an optimal 
and accurate exploitation of the next generation of galaxy surveys. 
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